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We study the topology and geometry of two dimensional coarsening foams with arbitrary liquid 
fraction. To interpolate between the dry limit described by von Neumann's law, and the wet 
limit described by Marqusee equation, the relevant bubble characteristics are the Plateau border 
radius and a new variable, the effective number of sides. We propose an equation for the individual 
bubble growth rate as the weighted sum of the growth through bubble-bubble interfaces and through 
bubble-Plateau borders interfaces. The resulting prediction is successfully tested, without adjustable 
parameter, using extensive bidimensional Potts model simulations. Simulations also show that a self- 
similar growth regime is observed at any liquid fraction and determine how the average size growth 
exponent, side number distribution and relative size distribution interpolate between the extreme 
limits. Applications include concentrated emulsions, grains in polycrystals and other domains with 
coarsening driven by curvature. 



Liquid foams, namely gas bubbles separated by a con- 
tinuous liquid phase, are ubiquitous [HE]. In floating 
foams as beer heads, ocean froths or pollutant foams, 
the fraction (j) of their volume occupied by the liquid de- 
creases with height, varying from a dry foam at the top 
to a bubbly liquid at the foam-liquid interface. 

Since pressure can differ from one bubble to another, 
gas slowly diffuses. Some bubbles disappear and, as no 
new one is created, the average size increases. Foam 
coarsening is analogous to that of concentrated emul- 
sions, grains in polycrystals, or two-phase domains where 
interface dynamics is driven by curvature. Its dynamics 
depends mainly on ^, up to a material-specific time scale 
determined by the foam physico-chemistry [1, 2 . 

Understanding foam coarsening requires two different 
levels. First, the individual bubble growth law^ which rules 
a bubble's growth rate according to its size or shape. This 
law can be stated as a static geometry problem, and may 
be obtained analytically or by detailed bubble shape sim- 
ulation. Second, the effect of such individual growth on 
the statistics of the foam^ i.e., bubble size and topology 
distributions, requires statistical theories or large bubble 
number simulations. 

In the very dry limit ^ ^ 0, bubbles are polyhedra with 
thin curved faces meeting by three along thin lines called 
Plateau borders. Coarsening in that limit has been inves- 
tigated experimentally, numerically and theoretically in 
two (2D) [3Hz] and later in three dimensions (3D) [8Ul4]. 

In the very wet limit ^ ^ 1, bubbles are round, dis- 
persed in the liquid and far from each other, forming a 
"bubbly liquid" rather than a foam stricto sensu. Their 
coarsening follows "Ostwald-Lifschitz-Slyozov- Wagner" 



ripening, in 3D [151 [16] and later in 2D fT7l[T8]. 

In both limits, the foam eventually reaches a self- 
similar growth regime: statistical distributions of face 
numbers and relative sizes become invariant. Only the 
average size (R) grows in time, as a power law (R) ~ t^, 
with /3 = 1/2 in the dry limit and /3 = 1/3 in the wet 
one, reflecting that the underlying physical processes are 
different. The number N of bubbles thus decreases as 
^-2/5 ^^ 2D and t~^^ in 3D. The growth law for interme- 
diary liquid fractions has been addressed in experiments 
[19 and simulations [20l[2T], but still lacks a unified the- 
oretical description. 

Here we address the 2D case. We propose a growth law 
to interpolate for < ^ < 1, with two parameters (dif- 
fusion coefficients) which are determined in each limit. 
To test our prediction on 2D foam coarsening experi- 
ments is difficult, because we are not aware of any study 
where is systematically varied and precisely measured, 
or even rigorously defined. We rather use numerical sim- 
ulations based on Potts model, suitable for large bubble 
numbers [10 11221124] . Beside testing our prediction, sim- 
ulations also show that, for any (/), side number and rela- 
tive size distributions reach a self-similar growth regime 
where (R) grows as t^ . Values of (3 interpolate between 
1/2 and 1/3. 

In the 2D dry limit, gas diffuses through neighbor bub- 
bles walls, due to the pressure difference between bubbles 
related with wall curvature. The walls are curved because 
they meet at threefold vertices, forming equal angles of 
27r/3. Each vertex is responsible for a turn of 7r/3 in the 
vector tangent to the bubble perimeter. Consequently 
bubbles with 5 sides or less are convex, while bubbles 



with 7 sides or more are concave, so that the wahs curva- 
ture plus a turn of 7r/3 at each vertex sums up to 27r [1] . 
The resulting growth dynamics is von Neumann's law [4]: 
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where a* and n* are, respectively, area and number of 
sides of the i^^ bubble; t is time; D^ depends on the 
foam composition and is expressed in m^s~^ as a diffusion 
coefficient. Remarkably, the rhs of eq. ([l]) involves only 
the bubble's number of sides and not its size or shape. At 
any time, bubbles with n* < 6 shrink while bubbles with 
n^ > 6 grow. Since for topological reasons the average 
bubble number of sides is 6 [U |26. , eq. ^ is compatible 
with gas volume conservation in the whole foam. 

In the 2D wet limit, gas bubbles are dispersed in a 
liquid matrix. Dynamics is a consequence of pressure 
difference in the gas contained in a bubble or dissolved 
in the liquid. This pressure difference is proportional 
to the wall curvature, which for a circular bubble is the 
inverse of its radius, R^. Marqusee [17] wrote the growth 
law using only R^ : 



^ = 2,«.f: = 2,D„/m, 



wherel^^y is another diffusion coefficient-like constant as 
above, and 
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where KjS are j order modified Bessel functions of sec- 
ond kind; ^ is the screening length (roughly, the typical 
distance beyond which bubbles do not feel the influence 
of each other); and Rc is the critical radius for which 
there is no growth, calculated by imposing total gas vol- 
ume conservation, i.e. ^ XI ^ <^* = 0- At any given time, 
bubbles with radius smaller than Rc lose gas while those 
with radius larger than Rc gain gas. 

Interpolating between eqs. ^ and Q seems difficult 
because they use very different variables: number of sides 
n* and radius R^. However, both equations involve the 
product of curvature times length {i.e. angle. Fig. fl]) of 
the interfaces through which the gas diffuses. We pro- 
pose that for any (j) the bubble growth rate is simply the 
superposition of growth through the interfaces shared ei- 
ther directly with other bubbles or with Plateau borders. 
It can thus be calculated as the weighted average of eqs. 
(IT]) and ([2]) . The weights are fractions of 27r angle carried 
by dry or wet parts of the i^^ bubble perimeter, 6J^ or 6^, 
which are sums of angles 6d or 6^ carried, respectively, 
by all dry or wet interfaces of the bubble (Fig. [l]), such 
that 0^ + Ql^ = 27r. We now detail how to perform this 
linear superposition, which turns out to be unexpectedly 
successful at all 6s. 




FIG. 1: (Color online) Angles Od or 0-w measure the change 
in direction carried by a dry or wet interface. Here Ueff is 
slightly larger than 5, the bubble is rather dry, Gd > Sw, 
Ru, < R. Each white region represents a Plateau border and 
is simulated as several liquid drops (not drawn here). 



In dry foams, at each vertex 6d = 7r/3, so that 6^ = 
nV/3. In wet foams, ©^^ = 27r. To interpolate, we char- 
acterize the bubble i by introducing its effective number 
of sides defined for any (j) as: 
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Although nlrr conveys no more information than Q^ or 

^ for a polygonal dry 



O^, its value is intuitive: 
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f2) bubble; n*^ <> = 6 for a circular wet bubble; 1 — n*^ p/6 or 



W^rr/Q are exactly the fractions of 27r carried by either 
dry or wet interfaces, respectively. Similarly, we char- 
acterize the bubble i by the curvature radius R^j of its 
Plateau borders; i?^ = for a dry bubble and i?^ = R^ 
for a wet bubble. With these variables n^rr^ R\j^ we can 
weight eqs. ([l]) and ([2| and interpolate for any (j): 

^ = \ [DdKff - 6) + D»<;//(i? j] . (5) 

To test eq. ([s]) we have implemented extensive Potts 
model simulations in the spirit of [10]. Like in experi- 
mental pictures, a simulation represents gas bubbles and 
liquid phase as connected regions on a square lattice of 
pixels. To ensure that liquid is evenly distributed be- 
tween Plateau borders (fig. IT]), and that (j) is conserved, 
we represent the liquid phase as a fixed number of tiny 
"liquid drops" of fixed area, and each Plateau border con- 
tains several of these drops. Each interface between gas 
bubbles represents a thin film made of two gas-liquid in- 
terfaces which mildly repel each other. We then assign to 
the interfacial energy between gas bubbles a value (2 — e) 
times the interfacial energy between gas bubble and liq- 
uid drop. Experimentally, in soap foams e is typically 
of order 10~^. Its precise value is not crucial in what 
follows, as long as e <C 1 (and £ > to avoid numerical 
unstabilities). Here we choose e = 0.01 as a compromise 
between realism and computing speed. Since liquid drops 
are free to cluster together, we assign to drop-drop inter- 
faces an energy 10^ times smaller. The liquid is sucked 
in the Plateau borders (for reviews see p!|l25]). 

There are JV = 8944^ sites; each site s is assigned with 
a label, S{s). There are N bubbles (5* = 1 to N) and 



A^^ drops {S = N -\-l to N -{- N^), assigned with a type 
r = 1 or 0, respectively. A configuration has an energy: 
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where s' stands for the sum over the first 36 neighbors of 
the site s, to avoid pinning to the grid [27 and to extend 
the range of the disjoining pressure up to three pixels; S 
and S' are the labels of sites s and s' ^ respectively; S is 
the Kronecker symbol; r and r^ are the types of 5* and S' ^ 
respectively; J(l, 1) = 1.99 J(1,0), J(1,0) = J(0, 1) = 
0.7, and J(0, 0) = 0.001 are the interfacial energies; a{S) 
is the current area of liquid drop S; at = 8 pixels is a 
target area common to all drops; A = 9 penalises any 
deviation from at. 

Simulations begin with A/'o = 2 x 10^ gas bubbles ran- 
domly dispersed over the grid, with smooth interfaces 
and a normal distribution of areas around the average 
^. For ^ > 0, in the initial configuration each vertex 
contains at least one drop. The total number of drops 
and the initial average area of gas bubbles are set ac- 
cordingly to the desired (j). 

The simulation dynamics follows Monte Carlo method. 
We randomly choose a site, temporarily change its label 
to the value of one of its neighbors, and calculate the 
change in energy AE. If A^ < this relabeling is ac- 
cepted. If AE > the change is accepted with probabil- 
ity exp(— A£^/T), where T is the fluctuation allowance, 
here taken as T = 3 to escape possible metastable states. 

We measure diffusion coefficients as follows. We first 
perform a simulation in the dry limit (/) = 0. Plotting 
the rhs and Ihs of eq. (IT]) determines by linear regres- 
sion Dd = 1.68. We then perform a simulation in the 
wet limit (j) = 0.9. Plotting the rhs and Ihs of eq. ([2| 
determines by linear regression D^j = 0.78. With these 
parameters, together with area and wet and dry perime- 
ters for each bubble from simulations, eq. ([5| predicts 
without any adjustable parameter each bubble growth 
for any intermediary (f). 

We measure angles on a square grid as follows. The 
simulation is halted. We add over all dry interfacial pixels 
of bubble i the probability that it would grow (or shrink) 
over (yielding to) other bubbles, calculated by the Monte 
Carlo method. This determines the growth rate of bubble 
i through its dry interfaces, GJ^, and thus GJ^ = D^G^^. 

pi 

We then obtain 0J^ = 27r — 0^, and i?^ = -^ where 

P^ is the wet interface length. During this measure, no 
change is performed. The simulation then resumes. 

We run simulations with ^ = 0, 0.02, 0.06, 0.18, 0.36, 
0.54, 0.72 and 0.90. For any 0, the evolution of the gas 
bubble number has a power law behavior (Fig. [2|a)), 
which is compatible with self-similar growth regimes. 
Fig. [2jb,c) shows that the power law exponent P varies 
as P{(j)) ~ 1/2 — (/)^'^/6. Thus P{(j)) decreases continuously 
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FIG. 2: (Color online) (a) Evolution of the number of gas 
bubbles for different values of 0, in log- log plot. (b,c) Power 
law exponent ^ versus 0, in (b) linear and (c) log-log scales. 

The red line is /3 = 1/2 - ^^-^^V^- 
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FIG. 3: (Color online) Snapshots after 20,000 MCS for = 0, 
0.02, 0.06, 0.18, 0.54, and 0.90. 



from 1/2 to 1/3, the expected limit values, with d/3/d0 
diverging at (/) ^ 0. It would be interesting to explain 
theoretically this variation of p with (j). 

Fig. [3] shows snapshots for different (j)s after 20,000 
Monte Carlo Steps (MCS). The liquid accumulates at 
the vertices for small 0s, and, as it increases, liquid 
goes also between bubbles. The ratio c^/^, where d 
the typical distance between bubbles, estimated as (i = 
2{^J\f /{Nit) — Y^(a)/7r), and ^ the screening length (eq. 
|2|, increases with (Fig. SI of [28 ). At = 0.54, 
d/^ > 1, and at = 0.90 the bubbles are not touching 
each other. 

Fig. [4] presents the average area growth rate of gas 
bubbles versus R^^ and in insets the distribution func- 
tion of Ryj/{Ryj). The superposition of plots taken at dif- 
ferent times indicates a self-similar growth regime. The 
agreement with the theoretical prediction (eq. [5| is ex- 
cellent. For = all bubbles have R^ = and this plot 
does not convey any information. For (j) = 0.02 and 0.06 
the noise arises from measuring the curvature radius of 
Plateau borders. For different values of i?^, distributions 
of bubble growth rates are presented in Fig. S2 of [28] . 

Fig. [5] presents the average area growth rate of gas 




FIG. 4: (Color online) Gas bubbles growth rate versus 
Rw/{Rw) at different times of the self-similar growth regime 
for = 0.02 to 0.90. Insets: distribution function of 

Rw/{Rw)- 
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FIG. 5: (Color online) Gas bubble growth rate versus 
Tieff/{neff) at different times of the self-similar growth 
regime for == to 0.54. Insets: distribution function of 

^eff/{neff). 



bubbles versus the effective number of sides, rieff. Again, 
the agreement with the theoretical prediction (eq. [5| is 
excellent. For (j) = 0.90 all bubbles have ng// = 6 and 
this plot does not convey any information. For different 
values of ng//, distribution functions of bubble growth 
rates are plotted in Fig. S3 of ^28j. Plots of growth rates 
versus neff/{neff) (Fig. S4 of [28]) and versus Ry^/ {Ry^) 
(Fig. S5 of [28 ) discriminate the relative contributions 
of dry or wet interfaces to the growth. 
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FIG. 1: Ratio of average distance between air bubbles and screening length ^ for different liquid fractions during scaling regime. 
Observe that, for (f) > 0.54, d > ^. The red line is a linear fit. 
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FIG. 2: Normalized bubble growth probability distributions, for the liquid fractions presented in the paper (0.0, 0.02, 0.06, 
0.18, 0.54, and 0.9), in the three distinct ranges of Rw/{Rw) shown in the upper left corner. Notice that, for = 0.0, (Rw) — 0. 
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FIG. 3: Normalized bubble growth probability distributions, for the liquid fractions presented in the paper (0.0, 0.02, 0.06, 
0.18, 0.54, and 0.9), in the ranges of ne///(ne//) presented in the lower right corner. Notice that, as increases, the distributions 
narrow around rie//, so no bubbles are found far from this rieff value. For = 0.54 and = 0.9 there are only bubbles with 
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FIG. 4: Area growth rate of air bubbles through the dry interfaces as a function of ^e///(^e//) for different Uquid fractions 
and at different instants of the scaUng regime. The insets present the probabihty density of ^e///(^e//) at the same instants. 
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FIG. 5: Area growth rate of air bubbles through the wet interfaces as a function of Rw/{Rw) for different Uquid fractions and 
at different instants of the scahng regime. The insets present the probabihty density of Rw/{Rw) at the same instants. 



